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Spin relaxometry of single nitrogen- vacancy defects in diamond nanocrystals 
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We report an experimental study of the longitudinal relaxation time (Ti) of the electron spin 
associated with single nitrogen- vacancy (NV) defects hosted in nanodiamonds (ND). We first show 
that Ti decreases over three orders of magnitude when the ND size is reduced from 100 to 10 nm 
owing to the interaction of the NV electron spin with a bath of paramagnetic centers lying on the 
ND surface. We next tune the magnetic environment by decorating the ND surface with Gd^^ 
ions and observe an efficient Ti -quenching, which demonstrates magnetic noise sensing with a single 
electron spin. We estimate a sensitivity down to ^ 14 electron spins detected within 10 s, using a 
single NV defect hosted in a 10-nm-size ND. These results pave the way towards Ti -based nanoscale 
imaging of the spin density in biological samples. 

PACS numbers: 



The ability to detect spins is the cornerstone of mag- 
netic resonance imaging (MRI), which is currently one 
of the most important tools in life science. However, 
the sensitivity of conventional MRI techniques is limited 
to large spin ensembles, which in turn restricts the spa- 
tial resolution at the micrometer scale (TJ [2]. Extend- 
ing MRI techniques at the nanoscale can be achieved at 
sub-Kelvin temperature with magnetic resonance force 
microscopy, through the detection of weak magnetic 
forces [T, T . Another strategy consists in directly sens- 
ing the magnetic field created by spin magnetic moments 
with a nanoscale magnetometer. In that context, the 
electron spin associated with a nitrogen-vacancy (NV) 
defect in diamond has been recently proposed as an ul- 
trasensitive and atomic-sized magnetic field sensor [5 . In 
the last years, many schemes based on dynamical decou- 
pling pulse sequences have been devised for sensing ac 
or randomly fluctuating magnetic fields with a single NV 
spin i6j-t9j . These protocols recently enabled nuclear mag- 
netic resonance measurements on a few cubic nanometers 
sample volume [TOl [11] and the detection of a single elec- 
tron spin under ambient conditions \TT . 

An alternative approach for sensing randomly fluctu- 
ating magnetic fields - i.e. magnetic noise - is based on 
the measurement of the longitudinal spin relaxation time 
(Ti) of the NV defect electron spin. Using an ensemble 
of NV defects and a Ti -based sensing scheme, Steinert et 
al. recently demonstrated magnetic noise sensing with a 
sensitivity down to 1000 statistically polarized electron 
spins, as well as imaging of spin-labeled cellular struc- 
tures with a diffraction-limited spatial resolution (^ 500 
nm) [13 . Bringing the spatial resolution down to few 
nanometers could be achieved by using a single NV defect 



integrated in a scanning device, e.g. with a nanodiamond 
(ND) attached to the tip of an atomic force microscope 
( AFM) [m [15] . With this application in mind, we study 
here the Ti time of single NV defects hosted in NDs, as a 
function of ND size and magnetic environment. We first 
report a decrease of Ti over three orders of magnitude 
when the ND size is reduced from 100 to 10 nm. This 
behavior is explained by considering the interaction of 
the NV spin with a bath of intrinsic paramagnetic cen- 
ters lying on the ND surface. We next tune the magnetic 
environment by decorating the ND surface with param- 
agnetic molecules. As expected, a strong Ti-quenching 
is observed when the surface spin density is increased. 
From our data, we estimate a sensitivity of Ti -based re- 
laxometry down to ~ 14 electron spins detected within 
10 s, using a single NV defect hosted in a 10-nm ND. 

The NV defect ground state is a spin triplet (5* = 1) 
with a zero- field splitting D = 2.87 GHz between a sin- 
glet state rris = and a doublet rris = ±1 [Fig. [lla)]. 
Owing to spin-dependent intersystem-crossing (ISC) to- 
wards intermediate singlet states, optical pumping leads 
to an efficient spin polarization into the rus = spin 
sublevel, while the spin state can be readout through 
spin-dependent photoluminescence (PL) [16 . These two 
properties enable the measurement of the Ti relaxation 
time of the NV defect electron spin by using the simple se- 
quence depicted in the inset of Fig.jljb). After initializa- 
tion into the rris = spin sublevel with an optical pulse, 
the NV defect is kept in the dark for a time r, causing the 
system to relax towards a mixture of states rUg = 0, ±1. 
The resulting electron spin state is readout by applying a 
second optical pulse. For a sufficiently short integration 
time (300 ns in this work), the readout PL signal X{r) 
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FIG. 1: (a) -Simplified energy-level structure of the NV defect. 
The zoom indicates the energy levels and transition rates used 
for studying the NV defect spin dynamics in the dark, (b)- 
Integrated PL signal as a function of the dark time measured 
for a single NV center hosted in a 20-nm ND. The solid line is 
a fit to Eq. ([l]), which yields Ti = IGibl /is and Tm = 160 ±30 
ns. The experimental sequence used to measure T\ is shown 
in inset. Laser pulses with 3-/is duration are used both for 
initialization of the NV defect in nis — and for spin-state 
read-out by recording the PL signal in a detection window 
corresponding to the first 300 ns of the optical pulses. 



can be written as X(r) ^ Aono(r) + Ai[n+i(r) + n_i(r)], 
where Aq and A\ < Aq are the PL rates associated 
with spin states m^ = and rris = ±1, respectively, 
and no,±i(r) are the spin populations before applying 
the readout optical pulse. These populations are eval- 
uated within the simplified four-level model shown in 
Fig. [ija), which includes the ground state spin sublevels 
rris = 0, ±1 and the lowest-lying singlet state, thereafter 
referred to as the metastable state. We define Ti as the 
decay time of the population no, hence 1/Ti = 3/coi, 
where /cqi is the two-way transition rate between rus = 
and m^ = ±1. At short time scale, the spin populations 
are also affected by relaxation from the metastable state 
which decays towards the ground state spin sublevels as 
nm{r) = n^(0)e"^/^'-, where T^ = {kmo + 2kmi)~^ is 
the metastable state decay time. The value of this param- 
eter is ^ 200 ns [16l [17] . Using classical rate equations 
within this four-level model, the PL signal X(r) can be 
written [TSl 



I{r) = X(oo) 
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The expressions of X(oc), Cm and Ci are given in the 
Supplemental Material [18]. A typical measurement of 
X(r) is shown in Fig. [lib) for a single NV defect hosted 
in a 20-nm-size ND, together with a fit to Eq. (IT]) yield- 
ing Ti = 16 ± 1 /is. This value is almost two orders of 



magnitude smaller than the one measured for single NV 
defects hosted in bulk diamond samples [T9| . 

To understand this behavior, the Ti time was stud- 
ied as a function of the ND size. We started from 
commercially available NDs (SYP 0.05 and 0.25, Van 
Moppes SA) produced by milling type-lb high-pressure 
high-temperature (HPHT) diamond crystals with a high 
nitrogen content ( [N] ^ 200 ppm) . The formation of N V 
defects was carried out using high energy (13.6 MeV) 
electron irradiation followed by annealing at 800° C under 
vacuum. The irradiated NDs were then oxidized in air 
at 550° C during two hours in order to remove graphitic- 
related defects on the surface and produce stable NV 
defects [20 . The NDs were finally spin cast on a glass 
cover slip and studied using a scanning confocal micro- 
scope combined with an AFM (Attocube Systems), all 
operating under ambient conditions [18 . For each pho- 
toluminescent ND, the PL intensity autocorrelation func- 
tion was first recorded in order to verify that a single NV 
defect was hosted by the crystal. For a set of single NV 
defects in isolated NDs, the Ti time was measured by fit- 
ting the relaxation curve X{t) to Eq. (IT]) and AFM mea- 
surements were used to infer the ND diameter do, defined 
as the maximum height in the AFM scan [Fig.l2]^a)]. The 
relaxation rate 1/Ti is plotted as a function of the ND 
size in Figure WCo) for a set of 51 single NV defects in 
isolated NDs with d^ ranging from 7 nm to 88 nm. An 
increase of the relaxation rate over three orders of magni- 
tude is observed when the ND size decreases. Indeed, Ti 
ranges from a few /is for the smallest (< 10 nm) NDs to 
up to 1 ms for the biggest ones (> 60 nm) [Fig.l2]^a)&(b)]. 
For NV defects hosted in bulk diamond samples, 
phonon-assisted processes are the main causes of longi- 
tudinal spin relaxation at room temperature, with Ti ly- 
ing in the 1-10 ms range [19 . Relaxation induced by 
paramagnetic impurities like nitrogen atoms (PI cen- 
ters), which are the most abundant paramagnetic defects 
in type-lb diamond, dominates only at low temperature, 
and results in Ti times that can be as long as 100 s at 4 K 
[19 . In NDs, a bath of paramagnetic centers covering the 
surface provides an additional channel for Ti relaxation. 
These impurities have been identified by numerous stud- 
ies [2TH26] and are mainly ascribed to dangling bonds 
with unpaired electron spins. For NDs with an oxygen- 
terminated surface, as those used in this work, Tisler 
et al. determined a density of surface spins a ?^ 1 — 10 
spin/nm^ using spin coherence measurements and ensem- 
ble EPR measurements [25 . We attribute the shortening 
of Ti of single NV spins in NDs to these surface paramag- 
netic centers (SPCs), which adds a contribution k^^^ to 
the transition rate k^^^^ of the bulk material, such that 
the overall rate is /cqi = k^^^^ + ^01^ • 

This hypothesis is tested by modeling the ND as a 
sphere and the SPCs as an ensemble of randomly fluctu- 
ating spins with a surface density a [Fig.[2]^c)]. The SPCs 
produce a fluctuating magnetic field B(t) with zero- mean 
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FIG. 2: (a)-Relaxation curves measured for a single NV defect hosted in a small ND (do = 13 it 3 nm, blue curve) and a larger 
one (do = 73 ± 10 nm, red curve). Data fitting with Eq. (IT]) (solid lines) gives Ti values of 3.6 ± 0.3 /is and 802 ± 136 /is, 
respectively. The corresponding AFM images are shown on top of the graph. (b)-Longitudinal spin relaxation rate 1/Ti of the 
NV defect electron spin as a function of the ND diameter. We note that no obvious correlations with the size was observed 
for the other fitting parameters in Eq. ([l]), which were found to be Cm = 0.084 ± 0.037, Ci = 0.21 ± 0.12 and Tm = 198 ± 72 
ns (mean ± s.d.). (c)-The ND is modeled as a sphere with diameter do with a bath of randomly fluctuating surface spins with 
density a. In (b), the markers are experimental data while the lines are the results of the calculation using this model for a NV 
spin located at the center of the sphere (solid line) and 3 nm below the surface (dotted line). The parameters of the calculation 



are T^ 



2 ms and a = 1 nm 



m 



{B{t)) =0, that is characterized by the spectral densities 
Sb^^) = I-^ Bk{t)Bk{t + r)e-*^^dr, where the three 
components k = x^y^z are assumed to be uncorrelated. 
For a central NV spin S = 1 with an intrinsic quantiza- 
tion axis along z, one has [27] 



JLSPC 



7e 



[SbS^q) + 5'b^(cjo)] 



(2) 



where 7e is the electron gyromagnetic ratio and ujq = 
27rD is the electron spin resonance (ESR) frequency of 
the NV defect. As highlighted by Eq. Q, longitudinal 
spin relaxation is caused by the transverse components 
of the magnetic noise at the ESR frequency of the cen- 
tral spin. Assuming correlation functions of the form 
{Bk{0)Bk{r)) = (5^)e~l^l/^^ where Tc is the correlation 



time of the magnetic field and {B'^) 
laxation rate reads 



1 



1 



o-ibulk 



HBl 



its variance, the re- 



(3) 



1 + ^nr. 



2^2 



Here Bj^ = {B'^) + (By) is the variance of the transverse 
magnetic field, and l/T^""^^ = 3k}^f^. Assuming S' = 1/2 
surface spins, the variance 5^ is calculated by summing 
the dipolar field at the NV's location from each randomly 
oriented SPC and the correlation time Tc is evaluated by 
considering intra-bath dipolar coupling ^18|. Since B^ 
depends on the exact location and orientation of the NV 
defect inside the ND, the calculation is performed for two 
extreme configurations [Fig. [2]^c)]. In the best-case sce- 
nario, the NV defect is located at the center of the sphere, 
while in the worst- case scenario, it is lying 3 nm below 
the surface - near the known photostability limit of the 



NV defect [28 - with its axis being parallel to the sur- 
face [18]. As shown in Fig. |2]^b), the results of the model 
capture fairly well the experimental data with <j = 1 
nm~^ [29 . More precisely the relaxation rate scales as 
1/c^o which stems from the l/d^ dependence of the spin- 
spin interaction integrated over a surface. This effect is 
responsible for the variation of Ti over several orders of 
magnitude when the size of the ND decreases [18] . 

In view of testing the ability of Ti relaxometry to 
detect changes in the local magnetic environment, the 
ND surface was decorated with additional paramagnetic 
species. This was achieved by spin casting an aqueous 
solution of gadolinium perchlorate molecules Gd(C104)3, 
containing paramagnetic Gd^+ ions (5* = 7/2), which is 
a well-known relaxation contrast agent in MRI. The Ti 
time was measured for a set of 33 single NV defects in iso- 
lated NDs (i) before any treatment, (ii) after a first treat- 
ment with 1 mM of Gd^+ solution, and (iii) after a second 
treatment with 10 mM. The substrate was patterned with 
a metallic grid for precise and repeatable identification of 
each individual ND over repeated treatment steps. The 
histograms of the measured 1/Ti rates are shown in Fig. 
IsJ^a). The distribution is clearly shifted towards higher 
relaxation rate after each treatment step, which indicates 
that single NV defects feel the magnetic noise induced by 
the external Gd^+ ions. In Fig. p]^b), the relaxation rate 
after the first and second treatment step is plotted as a 
function of the rate in the bare nanocrystal, i.e. before 
any treatment. Almost all the investigated NV defects 
undergo a significant decrease in their Ti . The quenching 
ratio T] = Ti^bare/^i, treated IS fouud to be 7^ ^ 7 ou aver- 
age after the first treatment with 1 mM of Gd^+ solution 
[Fig. |3]^b)]. From this value, we estimate a surface den- 
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FIG. 3: (a)-Histograms of the 1/Ti relaxation rate obtained 
from a set of 33 single NV defects hosted in isolated NDs. 
The measurement is performed before any treatment (top 
panel), after adding 1 mM of the Gd"^^ solution (middle 
panel) and after further adding 10 mM of solution (bottom 
panel). (b)-Relaxation rate measured after the first (circles) 
and second (triangles) treatment step as a function of the 
rate of the bare ND. The solid lines are data fitting with 
linear functions whose slope indicates the average quenching 
ratio T] = Ti,bare/Ti, treated- Wc obtalu T] ^ 7 (grccu line) and 
?7 ~ 31 (orange line). A dashed line of slope 77 = 1 is plotted 
for reference. 



sity of Gd^+ spins <jGd ^ 4 nm~^, corresponding to the 
detection of ^ 1000 spins for a 10-nm ND [18]. After the 
second treatment step with 10 mM of Gd^+ solution, we 
obtain on average 77 ^ 31 corresponding crcd ^ 70 nm~^. 
The dispersion in the measured values for r] is attributed 
to a non uniform surface spin density (Jq^. In particular, 
NV defects located close to the diamond-substrate inter- 
face should be less affected by the Gd^+ treatment. It 
is worth mentioning that a similar environment-induced 
quenching effect could be observed on the spin coherence 
time T2 as well, with however much smaller quenching ra- 
tios since T^""^^ <C T^^^^^ [13 [30]. Furthermore, Ti-based 
sensing schemes do not require coherent manipulation of 
the NV electron spin with microwave pulses. 

By applying a few more treatment steps with the 
Gd^+ solution, we then analysed the regime of strongly 
fluctuating magnetic environment, bringing Ti in the 
sub-microsecond range. As shown in Fig.|4] Ti-quenching 
is accompanied by a significant reduction of the Ti de- 
cay contrast, defined as Cf^ = max[X(r)]/X(oo) — 1. 
In the inset of Fig. [4| we plot Cf^ as a function of Ti 
together with the calculation based on a rate equation 
model that takes into account the full dynamics of the 
NV defect [181 EI] • The contrast reduction is mainly due 
to the overlap between Tm and Ti decays [see Eq. ([l])]. 
In addition, when Ti ~ T^, optical initialization in the 
rris = state and spin state readout become less efficient, 
thus reducing further the contrast. 

Based on our experimental results, we finally estimate 
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FIG. 4: X(t) relaxation curves measured after repeated treat- 
ments with the Gd^^ solution. Solid lines are fit to Eq. ([l]). 
Inset: effective Ti contrast Cf as a function of Ti. The solid 
line is the result of the calculation (with no fit parameter) us- 
ing a rate equation model with the parameters of NV centers 
in bulk diamond [18^ . 



the sensitivity of Ti relaxometry to small changes in the 
magnetic environment. For that purpose, we consider 
an optimized single-r measurement by fixing r ^ Ti/2, 
which converts a modification of spin relaxation into a 
change of the PL signal X(r) with optimal signal-to-noise 
ratio ^l3l HE] . Assuming a photon shot noise limited sig- 
nal, the smallest number of additional surface electronic 
spins 6Afmin that can be detected by a single NV defect 
located at the center of a ND with size do is given by 



SMm 



WAt 



dtfi^) , 



(4) 



where V includes both the finite contrast of the Ti re- 
laxation signal and the rate of detected photons, At is 
the integration time and f{a) is a slowly increasing func- 
tion of the intrinsic density of surface spins a fTS^. As 
expected, it is crucial to use NDs as small as possible in 
view of sensing magnetic noise from external spins. For 
a single NV defect hosted in a 10-nm ND with <j = 1 
nm~^ - i.e. corresponding to Ti = 6.3 /is according to 
the above model [Fig. [2|b)] -, a typical photon counting 
rate IZ = 10^ s~^ under cw optical illumination and a Ti 
contrast Ci = 0.2, we find (5A/'min = 14 spins within 10 s 
of integration. This result highlights that Ti relaxome- 
try with a single NV spin hosted in a ND is a promising 
resource to probe nanoscale magnetic field fluctuations 
with a sensitivity down to a few electron spins, within a 
time scale that is compatible with scanning probe tech- 
niques. Such probes might find important applications in 
life sciences, e.g. to image the spin density in biological 
samples with an unprecedented spatial resolution. 
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SUPPLEMENTARY INFORMATION 

Experimental setup 

All measurements were performed using a combined AFM/confocal microscopy setup operating under ambient 
conditions. A detailed description of the experimental setup can be found in Ref. [1 . Single NV defects hosted 
in diamond nanocrystals were optically excited with a laser operating at the wavelength A = 532 nm and their 
photoluminescence (PL) in the 650 — 800 nm range was detected using an avalanche photodiode. Laser pulses were 
produced with an acousto-optical modulator (MT200-A0.5-VIS) with a characteristic rising time of 10 ns. For all 
experiments, the optical pumping power was set at 1 mW, corresponding to the saturation power of the NV defect 
radiative transition. 



Derivation of the relaxation curve X(r) 



We first consider the NV defect ground state as a closed three-level system composed of the three spin sublevels 
rris = 0, ±1 with populations no,±i. The spin sublevels rus = and rris — ±1 are coupled by two-way transition rates 
of strength /cqi [see Fig. l-(a) of the main article]. After initialization into the vfis = spin sublevel with an optical 



pulse, the population dynamics are given by 



1 

«o(r) = 3 + 

n±i(r) = -± 



no(0) 



-r/Ti 



n+i(0)-n_i(0) 



-r/(3Ti) 



no(0) 



-r/Ti 



(5) 
(6) 



where Ti = [3fcoi] "'" and ni(0) is the initial population of spin sublevel rris = i. 
As indicated in the main text, the readout PL signal I{t) can be written as 

I{t) = Aonoir) + Ai[n+i{T) + n_i(T)] , 



(7) 



where Aq and Ai are the PL rates associated with spin states rris = and rris = ±1. Owing to spin-dependent PL of 
the NV defect, we consider Ai < Aq. Using Eq. ([5|, the signal X(t) can then be written as 



X(t) = J(oo) 



-de- 



-r/Ti 



(8) 



with 



X(oo) 



Aq + 2Ai 

3 

Aq-Ai 

Ao + 2Ai 



[3no(0) - 1] 



As expected, the readout contrast Ci is proportional to the difference of spin-dependent PL rate {Aq — Ai), and to 
the amount of initial spin polarization into the rris = spin sublevel. 

We now take into account the lowest-lying singlet state of the NV defect, thereafter referred to as the metastable 
state with population n^ [see Fig. l-(a) of the main article]. This population decays towards rris = with a rate 
kmo = Oikm and towards m^ = ±1 with a rate kmi = (1 — a)A:m/2. The metastable decay time is therefore given by 
Tm = k^ = {kmO) + 2kmi)~^ ' Typically a ~ 0.5 and k^ = Tm ~ 200 ns ^iS^. Within this four-level model, the 
populations after a time r read 



rimir) = nm{0)e ^/^- , 
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Using Eq. ([7|), the readout PL signal can then be written as 



X(t) = X(oo) 

corresponding to Eq. (1) of the main paper with 

Ao 



l-C^e-^/^^+Cie-^/^i 



X(oo) = 



■2Ai 



3 

Aq-Ai 

Ao + 2Ai 
nm(0) . 



[3no(0)-l + n™(0)] , 



(9) 

(10) 

(11) 
(12) 



Note that we have assumed a ^ 1/3 in order to obtain a simple formula for Ci. 

Although Ci does not depend explicitly on Ti, the initial populations no(0) and n^(0) as well as the PL rates Aq 
and Ai do, as soon as 1/Ti is not negligible compared to all other transition rates involved in the initialization process. 
Since Tm is the longest decay time of the NV center's dynamics besides Ti , this means that Ci is constant as long as 
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FIG. 5: (a)-Rate equation model used to simulate the X(r) relaxation curves. (b)-Simulated X(t) curves for various Ti times, 
as described in the text. (c)-Effective Ti contrast measured for various single NV defects as a function of Ti, which was tuned 
by adding Gd^^ ions on the nanodiamond surface. The solid line is obtained from the simulated X(t) curves, according to the 
full model depicted in (a). (d)-Measured X(r) curves corresponding to NV A3. 



Ti ^ Tm- However, when Ti gets closer to T^, the longitudinal spin relaxation prevents efficient spin initialization 
in m^ = through optical pumping, i.e. no(0) -^ 1/3, as well as efficient spin state readout, i.e. {Aq — Ai) vanishes. 
This explains why the contrast of Ti decay strongly decreases when Ti ~ Tm [see Fig. 4 of the main paper]. 

In order to gain further insight on this contrast reduction of Ti relaxation curves, the NV defect dynamics was 
modeled using the seven- level model shown in Fig. [sja). A detailed description of this model can be found in Ref. [3]. 
We used classical rate equations to calculate the time evolution of the populations. The optical pumping parameter 
P was varied in time in order to simulate the exact sequence employed in the experiment [see inset in Fig. 1(c) of 
the main article]. We used P = 1 when the laser is on and /3 = when the laser is off. The PL was numerically 
integrated over the first 300 ns of the readout optical pulse, as in the experiment. Apart from /cqi = l/3Ti that was 
taken as a variable, the values of the other transition rates were taken from Ref. ^, namely kr = 65 /is~^, /com = H 
/is~^, kim = 80 /is~^, kmo = 3 /is~^ and kmi = 1.3 /is~^. The dark time evolution r was varied in order to simulate 
the relaxation curves X{r) for different values of Ti, as depicted in Fig. Isfb). From these simulations, the effective 
contrast Cf^ = max [X(r)] /X(oo) — 1 was inferred as a function of Ti, leading to the solid line shown in the inset of 
Fig. 4 of the main article. 

In Fig. [sjc), this line is plotted again, together with the experimental data from a set of 6 single NV centers in 
distinct diamond nanocrystals, including the data shown in Fig. 4 of the main article (NV B6). The model captures 
correctly, on average, the observed behavior of contrast versus Ti. However, we note that the investigated NV defects 
all exhibit a different maximum contrast. This feature can be explained by variations of the transition rates for 
different NV defects in nanocrystals, which depend on the local strain, the presence of nearby defects, etc. In Fig. 
[sfd) we plotted the X{r) relaxation curves measured for NV A3, where we observe the extreme case of a nearly flat 
response, from which the Ti time could not be determined with a meaningful uncertainty. 



Derivation of Ti for a single NV defect interacting Avith a bath of surface spins 



Equation (3) of the main article gives the 1/Ti relaxation rate of the central NV spin placed in a fluctuating 
magnetic field B{t) with zero mean, characterized by a variance B^ (orthogonal to the NV defect quantization axis) 
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FIG. 6: (a) Notations for the derivation of the variance B^ of the transverse magnetic field at the NV's location. The symmetry 
axis of the NV center is denoted z, while the XYZ reference frame serves to define the spherical angles (0, cj)) that describe the 
position of the bath spins on the sphere (vector u^ for spin i). The two frames xyz and XY Z can be chosen to be different for 
mathematical convenience (if the NV spin is not at the center of the sphere) when performing the integration, (b) Notations 
for the derivation of the fluctuation rate i^dip of the spin bath caused by spin-spin interactions. 



and a correlation time Tn, i.e. 
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(13) 
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where uoq = 27rD is the electron spin resonance (ESR) frequency of the NV defect. 

In this section, we calculate the two quantities B^ and Tc for a bath of electronic spins distributed on the surface 
of the nanodiamond, modeled as a sphere of diameter do. The density of surface spins is denoted a. 



Variance of the transverse magnetic field 



The variance of the transverse magnetic field B^ felt by the central spin can be obtained by summing over the 
contributions of all bath spins S^ according to Bj^ = ^^B"^-. The dipolar field radiated by a spin S^ located at 
position Ti on the surface is given by 



^[S,-3(S,.u,)u,)] , 



(14) 



where u^ = rj/r, [see notations in Fig. |6ja)]. Tracing over a purely mixed state, described by a density matrix 
p = 2S+1 -**-2S+i where I2S+1 is the identity matrix of size 2S + 1, we obtain 



?2 _ /d2 \ , /d2 \ _r^^f„ru2 , d2 M _ />07e^^ ^ ^^ 2 + Ssin^ Qj 



Bii = {Bii) + {Bii) = Tr{p(i?^,, + B^^)} 



Ait 



(15) 



where Cs = 2S-\-i ^m=-s ^ ~ ^ ^^^ ^* ^^ ^^^ angle between r^ and the z-axis, which is the quantization 
axis of the central NV spin. For a bath of surface density a, whose surface is described in spherical coordinates by 
r = r(^, (j)) [Fig. |6|a)], the summation over all the bath spins gives 
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(16) 



For a spherical surface of diameter d^ and a central spin located exactly at the center of the sphere, we have r = (io/2 
and a(6>, (j)) = 6>, which yields 



Bl 



f^MTU 






ttCs 



4 ' 



(17) 



Considering S = 1/2 surface spins, we find B± = 26 mT nm^ x ^32-. For instance, B± = 260 /iT-rms for a 10-nm 
nanocrystal with a = 1 nm~^. 






If the central spin is offset by 5r along the z-direction, then r(^, (j)) = \/ {do/2)'^ — Sr'^ sin 6-\-Sr cos and a(6>, 



in Eq. (16). An offset along the transverse direction, e.g. along the x-axis, can be taken into account using the same 
formula for r{0^ (j)) but a modified angle a(^, (j)) such that sin^(a(6>, 0)) = ||u^ x u^^lp = cos^ + sin^ 6>sin^ 0, instead 
of sin^(Qf(6>, 0)) = ||u^ X u^lP = sin^ 6. The latter case is the one considered to obtain the red dashed line in Fig. 2(b) 
of the main paper, where the NV defect is located 3 nm below the nanocrystal surface with its quantization axis {z) 
parallel to the surface. 

Correlation time of the hath 

The total fluctuation rate R = I/tc of the bath spins can be decomposed into two main contributions R = i^dip+^vib 
where i^dip is due to intra-bath dipolar coupling while i^vib is caused by intrinsic vibrational spin relaxation. 

An estimate for i^^ip can be obtained by summing the dipolar interactions of a given spin S^ with all other spins 



Sj of the bath, according to hR^ip = \/^j-^i{Hf,) where Hij is the magnetic dipolar interaction 

Hij = ^1^ [S. • S, - 3(S. • u,,)(S, • u,,)] . (18) 



Here Vij is the radius- vector between the two spins and Uij — Yij/vij.^ as depicted in Figure p[b). For purely mixed spin 
states, the density matrix describing the two-spin system is p = (25+1)2 ^2S+i^^2S+i' The quantity {H'fj) = Ti{pH'^j} 
is then given by 



im = i'^) '^ ■ (!!•) 



For a bath of surface spins with a surface density a, we consider 

V- 1 /"+°° 27rr , 

Y.-^''^ T^dr, (20) 

where rmin is the minimum ahowed distance between two surface spins - e.g. the lattice constant for a crystal. The 
final expression for i?dip reads 
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For a bath of electronic spins with S = 1/2 and using rmin = 0.15 nm - i.e. the nearest neighbour distance in diamond 
- we find R^ip = 11 ns~^ nm xa^/^. Hence i^dip = H ns~-^ = 27r x 1.8 GHz by using a = 1 nm~^. 

Several EPR studies on various types of nanodiamonds have shown that the vibrational spin relaxation i^vib of the 
SPCs is around 1 ns~^ [4jl5j or less [6l[7]. These values are much smaller than i^dip, as calculated above for a typical 
surface spin density a = 1 nm~^. Therefore, we can neglect the vibrational contribution to the overall fluctuation 
rate, so that R ~ i^dip- 

Ti relaxation as a function of the nanodiamond size 



Replacing Eq. (17) in Eq. (13), i.e. for a NV spin located at the center of a ND of diameter do, the longitudinal 



spin relaxation rate writes 






where R^cr) = I/tc ~ Rd\p{,^) is given by Eq. (21). It can be seen that the second term in Eq. (22), which dominates 



in NDs, scales as ^q . This dependence, which simply stems from the d^ of the spin-spin interaction integrated over 
a surface, is responsible for the variation of Ti over several orders of magnitude when the size of the ND decreases. 
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FIG. 7: (a,b) Longitudinal spin relaxation rate 1/Ti of the NV defect electron spin as a function of the ND diameter. The 
markers are experimental data while the lines are the results of the calculation for a NV spin located at the center of the sphere 
(solid line) and 3 nm below the surface (dotted line), for several values of the surface spin density a. We assumed Ti^^^ = 2 
ms. We considered a vibrational contribution to the spin bath fluctuation rate i?vib = in (a), as in the main article, and 
Rvih = 50 ns~^ in (b). 



In Fig. JTJ^a), we plot again the 1/Ti vs. do data points, together with the calculation using the above formulas. 
The calculation is performed for the two extreme cases where the NV spin is located either at the center of the sphere 



[best case, given by Eq. (22)] or 3 nm away from the surface with the NV axis parallel to the surface {worst case). The 
result of the calculation is shown for three different surface densities a. The value a ~ 1 nm~^ is the one that best 'fits' 
the experimental results, if the criterion is that most data points must lie in between the two extreme case theoretical 
curves. Figure mb) is identical to Fig. iTfa) but here we assumed a non-zero value for the vibrational contribution to 
the correlation time, namely i?vib = 50 ns~^, which is a typical value for electron spins in paramagnetic complexes [8], 
but is probably an extreme case for our nanodiamond SPCs. Now the experimental results are best 'fitted' with a 



surface density a ^ 
reasonable values. 



2 nm . This shows that the outcome of the model is not very sensitive to the choice of i^vib for 



Estimation of the number of gadolinium spins 

In this section we estimate the number of gadolinium spins that were added to the ND surface from the experi- 
mentally measured quenching ratio r] = Ti^bare/^i, treated [see Fig. 3 of the main text]. Here Ti^treated (resp. Ti^bare) 
denotes the relaxation time of the NV defect after (resp. before) adding Gd^+ ions. Since these experiments were 
performed in the limit 77 ^ 1, one can consider that Ti ^treated is only due to a bath of S' = 7/2 Gd^+ spins with a den- 
sity (JGd, producing at the NV location a fluctuating magnetic field characterized by a variance B'^ and a correlation 
time r^. The spin relaxation rate after adding Gd^+ ions can therefore be written 
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(23) 



Neglecting l/T^^^ik ^^^ considering only the intra-bath dipolar contribution to r^, one obtains from Eqs. (17) and 
(I2TI) that 
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In the latter expression, we assumed a = 1 nm~^, i.e. 1/tc = 11 ns~^ and l/r^ ^ ujq as justified below. Within 
this framework, rj depends neither on the ND size nor on its exact shape, in agreement with the experimental data. 
By adding 1 mM of Gd^+ solution, an averaged quenching ratio r^ = 7 is obtained [see Fig. 3 of the main text]. 
From this value, we estimate a surface density of gadolinium spins crcd ^ 4cr ~ 4 nm~^. For a 10-nm size ND, this 
corresponds to the detection of ~ 1000 gadolinium spins. After adding 10 mM of solution, we find 7^ = 31 on average, 
corresponding to crcd ^ 70cr ^ 70 nm~^. 

We note that crcd = 4 nm~^ yields i?^-p ^ 500 ns~^, which is indeed much larger than the vibrational contribution 
(typically i?(^-^ ~ 50 ns~^ [8 ) and than cjq = 18 ns~^. This validates the assumptions made to obtain Eq. (25), i.e. 
that only the intra-bath dipolar term contributes to r^ and 1/r^ ^ ujq. 

Sensitivity estimation 

In this section we derive the sensitivity of Ti relaxometry to a change in the number Me of external electronic spins 
on the surface, hence to a change in the relaxation rate Fi = 1/Ti of the NV spin. We consider a single-r detection 
scheme [9 that consists in measuring the PL signal X{t) after a fixed dark time r ^ T^, as shown in Figure [s] The 
number of photons detected while repeating this sequence during a total acquisition time At reads 

7V(ri)^^^At[l + Cie-^^^] , (26) 

where IZ is the photon counting rate under cw optical illumination. Tint is the width of the integration window (300 
ns in our experiments), and Ci <C 1 is the contrast. An infinitesimal change 5Ti of the relaxation rate is converted 
into a modification SN^^^^qx of the number of detected photons given by 

(^A^signai ^ (5Fl7^TintAtCle-^^" . (27) 

Assuming a photon shot noise limited signal (5A/'noise = \/N{Ti) and Ci <C 1, the signal to noise ratio (SNR) reads 

SNR = ^-^^ = ST, \/^7i;;;7AiCie-r^^ . (28) 

The maximal signal-to-noise ratio SNR^ is obtained for r = Ti/2 and is given by 



SNR. = i^Cn/^^^ . (29) 

V^ V 2e ^ ^ 

For a net change JFi, the SNR is therefore improved when Fi decreases. However STi and Fi are not independent 
quantities, as we shall see below. 

In the following, we consider a single NV spin located at the center of a spherical diamond nanocrystal of diameter 



do. Using Eqs. (21)-(22) and neglecting 1/T{^^^^, the relaxation rate due to a bath of electronic spins S = 1/2 with 



surface density a can be expressed as 

with A = 5.75 x 10^ s~^nm^ and B = 2.64 nm~^. We seek to calculate the effect of an infinitesimal increase Sa of 
the surface spin density. The resulting modification of the relaxation rate 6Ti is given by 

STi _Sa n B 



^^ = di7TB' (^«) 



(31) 



Inserting Eqs. (30) and (31) into Eq. (29), we finally obtain 



SNR™ = S- X ^V 2e .i/4(. + 5)3/2 • (32) 



In terms of increase in the number of surface spins 5JVe = ScrTrdQ^ the SNR writes 



SNR™ = <5A4 X ^V 2e <xV4(, + 5)3/2- (33) 
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FIG. 8: (a)-Sequence used in the single-r detection scheme. In the derivation of the sensitivity, the repetition period is taken 
to be equal to the dark time r, i.e. we neglect the duration Tpuise of the laser pulse, which can be as short as 1 /is in practice. 
(b)-Typical X(t) curve calculated with a spin relaxation rate Fi = 10^ s~^ (blue curve). The red curve illustrates the effect of 
an increase of the relaxation rate by ^Fi = 0.2Fi. The single-r detection scheme enables to probe the PL difference for a given 
dark time r, which is denoted ^A^signai in the text and is related to the change ^Fi. 



Owing to the d^ ^ dependence, it is crucial to use a nanodiamond as small as possible in order to detect a net change 
6J\fe. The smallest number of spins that can be detected is obtained for a signal to noise ratio of one, and reads 
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where V 



27r 



A7^Ti, 



and f{(j) 



(7^/^(a+B)^ 



2g ct^^v^ J y^ J — a-\-3B ^^ ^ slowly fluctuating function of a. For a single NV spin located 

at the center of a 10-nm diameter nanocrystal with an intrinsic spin density a = 1 nm~^, a typical photon counting 
rate 71 = 10^ s~^ and a contrast Ci = 0.2, we find 5Me,min ^ 14 spins within 10 s of integration. We emphasize that 
although the above calculation assumes the 5Me additional spins to be distributed all over the nanocrystal's surface, 
we expect our conclusions, especially the distance dependence, to be true even if the 5Me spins are localized in a 
point-like spot, e.g. in a molecule that would be approached close to the nanocrystal's surface. 
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